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Abstract 

We have measured the cross section a(e + e~ — > 7r + 7r~7) at an energy W = = 
1.02 GeV with the KLOE detector at the electron-positron collider DA$NE. From 
the dependence of the cross section on the invariant mass of the two-pion system, 
we extract a(e + e~ — ► 7r + 7r~) for the mass range 0.35 < s < 0.95 GeV 2 . From this 
result, we calculate the pion form factor and the hadronic contribution to the muon 
anomaly, a M . 



Key words: Hadronic cross section, initial state radiation, pion form factor, muon 
anomaly 
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1 Hadronic cross section at DA$NE 



1.1 Motivation 

The recent precision measurement of the muon anomaly at the Brookhaven 
National Laboratory [1] has led to renewed interest in an accurate measure- 
ment of the cross section for e + e~ annihilation into hadrons. Contributions to 
the photon spectral functions due to quark loops are not calculable for low- 
hadronic-mass states by perturbative QCD at low energy. However, they can 
be obtained by connecting the imaginary part of the hadronic piece of the po- 
larization function by unitarity to the cross section for e + e~ — > hadrons [2,3]. 
A dispersion relation can thus be derived, giving the contribution to as an 
integral over the hadronic cross section multiplied by a kernel K(s), which 
behaves approximately like 1/s: 

1 °° 

< d = ^ / <J e+e -_ hadl (s)K(s)ds. (1) 

Ami 

The process e + e~^7r + 7r~ below 1 GeV accounts for 62% of the total hadronic 
contribution [4]. The most recent measurements of cr(e + e~^7r + 7r~) for values 
of y/s between 610 and 961 MeV come from the CMD-2 experiment at VEPP- 
2M where the quoted systematic error is 0.6% and the contribution of the 
statistical error on ajf d is ~0.7% [5,6]. These data, together with r and e + e~ 
data up to 3 GeV, have been used to produce a prediction for comparison with 
the BNL result [7]. There is however a rather strong disagreement between the 
aj} ad value obtained using r decay data after isospin-breaking corrections and 
e + e~^7r + 7r~ data. Moreover, the e + e~^7r + 7r~ based result disagrees by ~3<r 
with the direct measurement of a,,. 



1.2 Radiative Return 

Initial state radiation (ISR) is a convenient mechanism that allows one to 
study e + e~^ hadrons over the entire range from 2m n to W, the center of mass 
energy of the colliding beams. In this case, there are complications from final- 
state radiation (FSR). For a photon radiated prior to the annihilation of the 
e + e~ pair, the mass of the 7r + 7r~system is 1 m w + n - = \JW 2 — 2WE 1 . Instead, 
for a photon radiated by the final-state pions, the virtual photon coupling to 
the 7r + 7r~ pair has a mass W . By counting vertices, the relative probabilities 

1 Neglecting the small momentum of the (/). 
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of ISR and FSR are of the same order. This requires careful estimates of the 
two processes in order to be able to use the reaction e + e~^7r + 7r~7 to extract 
o-(e + e~^7r + 7r~). The Karlsruhe theory group has developed the EVA and 
PHOKHARA Monte Carlo programs which are fundamental to our analysis 
[8-13]. In particular, the PHOKHARA Monte Carlo simulation has been used 
to evaluate the contribution for the ISR process (via the radiation function 
H) in order to derive the hadronic cross section: 

^ tfoVr+Tr 7 = ^^(s^His^), (2) 



where s n = m 2 + _, which coincides with the invariant mass s of the intermedi- 
ate photon for the case of ISR radiation only. The equation above is correct at 
leading order if FSR emission can be neglected. The case of NLO terms, with 
the simultaneous emission of ISR and FSR photons, is discussed in Sec. 3.1. 

The present analysis is based on the observation of Ref. [8] that for small 
polar angle 6 1 of the radiated photon, the ISR process vastly dominates over 
the FSR process. In the following we restrict ourself to studying the reaction 
e + e~^7r + 7r - 7 with # 7 < 15° or # 7 > 165°. For small s^, the di-pion system 
recoiling against a small angle photon will result in one or both pions being 
lost at small angle as well. We are therefore limited to measuring (j(7r + 7r~) 
for v /s^ >550 MeV. In the future extension of this work we will be able to 
measure the cross section near threshold. This is very important, since there 
are no good measurements of a(7r + 7r _ ) at low masses, which weigh strongly 
in the estimate of a^ ad . 



2 Measurement of a(e + e — > ty + tt 7) 



2. 1 Signal selection 



The KLOE detector consists mainly of a high resolution drift chamber with 
transverse momentum resolution <j Pt /pt < 0.4% [14] and an electromagnetic 
calorimeter with energy resolution o~e/E = 5.7%/y / £'(GeV) [15]. In the cur- 
rent analysis, we have concentrated on events in which the pions are emitted 
at polar angles 6 n between 50° and 130°. The direction and energy of the 
photon is reconstructed from the pion tracks by closing the kinematics; ex- 
plicit photon detection is not required. As a consequence, a requirement on 
the di-pion production angle 6 wn smaller than 15° (or greater than 165°) is 
performed instead of a requirement on the photon angle # 7 . The acceptance 
regions are shown in Fig. 1, left. These specific acceptance requirements re- 
duce background contamination and the relative contribution of final-state 
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Fig. 1. Left: schematic view of the KLOE detector with the angular acceptance 
regions for pions (horizontally hatched area) and photons (vertically hatched area). 
The photon angle is evaluated from the two pion tracks. Right: 2-dimensional re- 
quirement in the plane of m tr k/MeV and s T /GeV 2 . 

radiation from the pions to very low levels [16]. It will be shown in the follow- 
ing that an efficient and nearly background free signal selection can be done 
without explicit photon tagging. 

The selection of e + e~ — > 7r + 7r~7 events is performed with the following steps: 

• Detection of two charged tracks connected to a vertex: Two charged tracks 
with polar angles between 50° and 130° connected to a vertex in the fidu- 
cial volume, R xy < 8 cm, \z\ < 7 cm, are required. Additional requirements 
on transverse momentum, > 160 MeV, and on longitudinal momentum, 
\p z \ > 90 MeV, reject spiralling tracks and ensure good reconstruction con- 
ditions. 

• Identification of pion tracks: Separation of pions from electrons is performed 
using a PID method based on approximate likelihood estimators. These 
estimators are based on the comparison of time-of-flight versus momentum 
and on the shape and energy deposition of the calorimeter clusters produced 
by charged tracks. The functions have been built using control samples 
of 7T + 7r~7r° and e + e~7 events in data, in order to obtain the calorimeter 
response for pions and electrons. An event is selected as signal if at least 
one of the two tracks is identified as a pion. In this way, the content of e + e~7 
events in the signal sample is drastically reduced, while the efficiency for 
retaining 7r + 7r~7 events is still very high (> 98%). 

• Requirement on the track mass: The track mass (m tr k) is a kinematic vari- 
able corresponding to the mass of the charged tracks under the hypothesis 
that the final state consists of two particles with the same mass and one 
photon. It is calculated from the reconstructed momenta of the 7r + and 7r~ 
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(p+, p-) and the center-of-mass energy W. Requiring a value larger than 
120 MeV rejects /i + /i~7 events, while in order to reject 7r + 7r~7r° events, an 
s^-dependent requirement is used (see Fig. 1, right). 
• Requirement on the di-pion angle 9 n7V : The aforementioned requirement on 
the di-pion angle 9 W7T < 15° or > 165° and 50° < 9 n < 130° is performed. 
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Fig. 2. Distribution of counts as a function of s w , in bins of 0.01 GeV 2 , after applying 



the acceptance and selection requirements. C = 141.4 pb 1 ; data from 2001. 



The data used for the analysis were taken from July to December 2001, 
yielding an integrated luminosity of C = 141.4 pb _1 . After the selection re- 
quirements mentioned above, we find 1.555 x 10 6 events, corresponding to 
~ 11000 events/pb . Fi gure 2 shows the distribution of the ir + ii 7 events in 
bins of 0.01 GeV 2 for s n . The p peak and the p — uj interference structure are 
clearly visible, even without unfolding the spectrum from the detector resolu- 
tion, demonstrating the excellent momentum resolution of the KLOE detector. 
To obtain the cross section for 0° < 9 W < 180° and 9 n7r < 15°, 9 n7T > 165° we 
subtract the residual background from this spectrum and divide by the selec- 
tion efficiency, acceptance, and integrated luminosity: 

Gfov+7r- 7 = AiVo bs - AA^Bkg 1 1 ^ 



The background subtraction, the evaluation of the selection efficiency and the 
acceptance, the measurement of the integrated luminosity, and the unfolding 
of the experimental resolution on s n (omitted from Eq. (3) for clarity) are 
discussed below. Detailed information on all the aspects of the analysis is 
available in [17]. 
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2.2 Background subtraction 



After applying the requirements on the fiducial volume, the likelihood, and 
m tr k, a residual background of e + e~7 , /i + /i~7, and 7r + 7r~7r° events remains. 
The population of signal and background events in the [s^m^] plane is illus- 
trated in Fig. 1, right. Background from e + e~7 and yU + /i _ 7 events is concen- 
trated at low values of m tr k ■ The amount of background in the signal region 
is obtained by fitting the m tr k spectra of the selected events (except for the 
m tr k requirement) in slices of s n . The m tr k spectra for signal and /z + /i~7 events 
are obtained from Monte Carlo simulation, while for e + e~7 events, m tr k is ob- 
tained directly from data, using a dedicated sample of 152 pb _1 . An example 
of such a fit to determine the background fraction for /i + /i~7 events is shown 
in Fig. 3, left. Background from 7r + 7r~7r° events appears at higher m tr k values 
and the missing mass m^ iss = (p^ — p + — p_) 2 , peaks at m 2 ^. The number of 
7T + 7r~7r° events in the signal region is obtained by fitting the m miss distribu- 
tion with the shapes obtained from the Monte Carlo simulation; an example is 
shown in Fig. 3, right. The shape of the background distribution is well repro- 




Fig. 3. Left: fit of the mass peak for muons. Right: ir mass fit. These fits are used 
to estimate the /i + ^~7 and ir + iT~ir backgrounds to the 7r + 7r~7 channel. Points are 
data, solid line is Monte Carlo simulation. Dashed line and hashed area represent 
the Monte Carlo evaluation of background (/i + /i~7 in the left plot and 7r + 7r~7r° in 
the right one) and tt + tt~^ contributions, respectively. 

duced by the Monte Carlo simulation, ensuring that systematic uncertainties 
are smaller than the fit errors, which are considered as systematic errors of 
the procedure. 

The contribution of all backgrounds to the observed signal is below 2% above 
0.5 GeV 2 , while it increases up to ~ 10% at s n = 0.35 GeV 2 . Other possible 
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sources of background for which the contributions have been evaluated are the 
process e + e~ — > e + e~7r + 7r~ with the electrons emitted along the beam pipe [18] 
and the NLO corrections to the FSR in the process e + e~ — > /i + /i~7 [19]. The 
systematic uncertainties associated with the background estimates for all these 
sources have been added in quadrature; the results are shown in Table 2. 



2.3 The selection efficiency es e i 

The selection efficiency is the product of the efficiencies associated with the 
trigger, the event reconstruction, the background filtering and the track mass 
requirement: in 

• Trigger efficiency: Events in the 7r + 7r~7 sample must satisfy the calorimeter 
trigger, i.e., there must be at least two trigger sectors with energy deposi- 
tion above threshold (for details on the KLOE trigger see [20]). The trig- 
ger also includes a cosmic-ray veto: events with energy deposition above a 
certain threshold in the outermost layer of the calorimeter are rejected on- 
line. While the calorimeter trigger itself is rather efficient for signal events 
(> 95%), the cosmic-ray veto rejects a significant fraction of 7r + 7r _ 7 events 
since such events mimic cosmic rays. The cosmic-ray veto inefficiency is 
on the level of few percent at small values of s n < 0.4GeV 2 but reaches 
up to 30% at = 0.95GeV 2 . The overall trigger efficiency, including the 
effect of the cosmic-ray veto, was evaluated from the probability for sin- 
gle pions to fire trigger sectors in 7r + 7r~7 events wherein part of the event 
could be ascertained to have satisfied the trigger alone. The fractional un- 
certainty associated with this procedure was estimated to be <^6trg (^7r) = 
[exp(0.43 — 4.95,,-fGeV 2 ]) + 0.08] (expressed in percent), and is dominated 
by the systematics of establishing the correct track-to-trigger sector associ- 
ation. 

• Background filter efficiency: During reconstruction, an offline filtering pro- 
cedure identifies and rejects background events as soon as they have been 
reconstructed in the calorimeter [21]. The efficiency of this filter has been 
studied using a dedicated sample of 7r + 7r~7 events that were rejected by the 
filter itself. Since the filtering procedure is very sensitive to the presence of 
accidental clusters in the electromagnetic calorimeter, the efficiency of this 
filter was parametrized as a function of the background conditions during 
data taking and averaged over time. The filter efficiency was found to be 
uniform in s w and 96.6% on average, with a flat systematic error of 0.6%. 

• Tracking efficiency: The tracking efficiency (96% and uniform in s n ) was 
evaluated using 7r + 7r~7r° and 7r + 7r~events identified by calorimeter infor- 
mation plus the presence of one fitted track. The single-track efficiency as 
a function of p n and 9 n was compared with Monte Carlo simulation; the 
difference, on the order of 0.3% and flat in momentum, was taken as the 
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systematic error of this procedure. 

• Vertex efficiency: The efficiency of the vertex finding algorithm has been 
evaluated via Monte Carlo simulation and checked with a sample of 7i + n^n° 
and 7r + 7r~7 events obtained from data. The absolute vertex efficiency at low 
energies is 91% and is increasing up to 97% at high values of s n . An uncer- 
tainty of 0.3%, uniform in s n , is taken as the contribution to the systematic 
error for this efficiency. 

• Pion identification: The efficiency for -n/e separation has been evaluated 
by selecting 7r + 7r~7 events on the basis of one track and examining the 
distribution of the likelihood estimator for the other one. In the analysis, 
only one track is required to satisfy the likelihood requirement, for which 
the efficiency is > 99.9%. Therefore, no correction for the inefficiency on 
pion identification needs to be applied; the contribution to the systematic 
error is taken to be 0.1%. 

• Track mass: The efficiency of the m trk requirement is obtained as a by- 
product of the background evaluation; the result of the fit provides the effi- 
ciency in each s n bin. However, this efficiency depends upon the treatment 
of multi-photon processes in the Monte Carlo simulation. The m tr k effi- 
ciency has been obtained with our reference Monte Carlo simulation, which 
uses PHOKHARA. To check the efficiency determination we have compared 
PHOKHARA with BABAYAGA [22], which is the generator used for the 
luminosity measurement. In the latter generator, ISR is treated using the 
parton-shower approach. The resulting value for the m tr k efficiency differs 
from that evaluated with PHOKHARA by 0.2%. Effects on the efficiency 
from the simultaneous emission of an ISR and a FSR photon are discussed 
in Sec. 3.1. 



2.4 Unfolding of the mass resolution 

To obtain da n7T ^/ds n as a function of the true value of s n , we unfold the mass 
resolution from the measured s n distribution. The measured value of s^obs is 
related to the true value via the resolution matrix G(s vrit ruc ~ s^obsls^truc), 
which has been obtained by a Monte Carlo simulation carefully tuned to re- 
produce the data. The resolution matrix is nearly diagonal, as can be seen in 
in Fig. 4, left. A comparison of the track-mass distribution for data and Monte 
Carlo events is shown Fig. 4, right. 

Unfolding of the spectrum is performed using GURU [23], an unfolding pro- 
gram based on the singular value decomposition (SVD). We found that the 
systematic error due to unfolding is dominated by the uncertainty on the 
value chosen to regularize the procedure itself. Table 3 shows the systematic 
uncertainty as function of s n , introduced into the 7r + 7r~7 spectrum due to the 
unfolding. These values are taken to be correlated errors, and translate into a 
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Fig. 4. Left: Smearing matrix representing the correlation between generated 
(svr,true) and reconstructed (s^obs) values for s n ; the high precision of the DC re- 
sults in an almost diagonal matrix. Right: track mass distribution expected from 
the Monte Carlo simulation compared to the experimental one. 

0.2% systematic uncertainty on a M . 2 

In addition, the unfolding procedure correlates the statistical errors in the 
7T + 7r _ 7 spectrum (see also Sec. 4). 



2.5 Acceptance correction 



After all corrections discussed above, we obtain the spectrum for 7r + 7r~7 events 
denned by the acceptance requirements 50° < 9 n < 130°, 9 nn < 15° or 
9 nn > 165°, pt > 160 MeV, and p z > 90 MeV. To derive the cross sec- 
tion for the process e + e~ — > 7r + 7r~7 with 9 n7T < 15° or 9 n7T > 165°, the effects 
of the other requirements on the momentum and polar angle of the pions 
have been evaluated using PHOKHARA. The systematic error of 0.3% on the 
acceptance fraction has been estimated by a comparison of data and Monte 
Carlo distributions. 



2.6 Luminosity measurement 



The integrated luminosity is measured with the KLOE detector itself using 
very-large- angle Bhabha (VLAB) events. The effective Bhabha cross section 
at large angles (55° < 9 + - < 125°) is about 430 nb. This cross section is 

2 This value should be considered as an overestimate of the real effect introduced 
by the unfolding procedure on a^, as discussed in [17]. 
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large enough so that the statistical error on the luminosity measurement is 
negligible. The number of VLAB candidates, -/Vvlab j is counted and normal- 
ized to the effective Bhabha cross section, <7vlab> obtained by Monte Carlo 
simulation, after subtraction of the background, <5Bkg : 



The precision of the luminosity measurement depends on the correct inclu- 
sion of higher-order terms in computing the Bhabha cross section. We use the 
Bhabha event generator BABAYAGA [22], which has been developed explic- 
itly for DA$NE. In BABAYAGA, QED radiative corrections are taken into 
account in the framework of the parton-shower method. The precision quoted 
is 0.5%. The result for the effective Bhabha cross section has been compared 
with that from BHAGENF [24,25], a full order-a event generator. We find 
agreement to better than 0.2%. 

VLAB events are selected with requirements on variables that are well repro- 
duced by the KLOE Monte Carlo simulation. The electron and positron polar 
angle reqiurements, 55° < 6 + _ < 125°, are based on the calorimeter clusters, 
while the energy requirements, E + _ > 400 MeV, are based on drift cham- 
ber information. The background from [i + fi~ (7), 7r + 7r~(7) and 7r + 7r~7r° events 
is well below 1% and is subtracted. All selection efficiencies (trigger, EmC 
cluster, DC tracking) are > 99% as obtained by Monte Carlo simulation and 
confirmed with data. We obtain excellent agreement between the experimental 
distributions (#+,_, E + _) and those obtained from Monte Carlo simulation, 
as seen in Fig. 5. Finally, corrections are applied on a run-by-run basis for 
fluctuations in the center-of-mass energy of the machine and in the detec- 
tor calibrations. The experimental uncertainty in the acceptance due to all 
these effects is 0.3%. We assign a total systematic error on the luminosity of 
5C = 0.5%th©0.3% CX p- The luminosity measurement is independently checked 
using e + e~— ^77 events. We find agreement to within 0.2%. 



2.7 -k + tt 7 cross section 



Our results for the differential cross section da(e + e — > 7r + 7r 7)/cis 7r with 



0° <6 n < 180° and 9^ < 15°, 6 W7T > 165° are plotted in Fig. 6, left, and are 
presented in numerical form in the second column of Table 1. 




(4) 
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Fig. 5. Data-Monte Carlo comparison of the (left) and E+- (right) distribu- 
tions for Bhabha events selected at large angle as described in the text. 



3 Extraction of a(e + e — > 7r + 7r ) and {F^s)] 2 



In order to extract the e + e~ — > 7r + 7r~ cross section, the radiation function H 
is needed (see Eq. [2]). This function is obtained from PHOKHARA, setting 
F n (s) = 1 and switching off the vacuum polarization of the intermediate pho- 
ton in the generator. Applying Eq. (2), and taking the FSR contribution into 
account, as described in the following section, the hadronic cross section as a 
function of the invariant mass of the virtual photon, s = m 2 t , is obtained, as 
shown in Fig. 6, right. 





Fig. 6. Left: differential cross section for the e + e — ► ir + ir 7 process, inclusive in 
9 W and with 9 n7r < 15° (9 W7T > 165°). Right: cross section for e + e~ — > 7r + 7r~. 
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3.1 FSR corrections 



Events with one or more photons emitted by the pions (FSR) without any 
photons in the initial state must be considered as a background to our mea- 
surement. Our event selection strongly suppresses the contribution of such 
events to well below 1% over the entire range of s n . 

However, events with the simultaneous emission of one photon from the initial 
state and one photon from the final state must be included in our selection in 
order for the e + e~ — > 7r + 7r~ cross section to be inclusive with respect to FSR 
(see Ref. [7] for details). More specifically, since the radiator function H only 
describes the ISR part of the radiative corrections, the process e + e~ — > 7* — > 
7r + 7r~7 ISR (7FSR), with one photon from initial state and possibly another from 
the final state, corresponds to e + e~ — > 7* — > 7r + 7r~(7 F sR.) after the division by 
H. 

Therefore, 

na 2 c fo- 7 ™ i "7(7FSR) 
„(eV - ^-(7 FSR )) = 17^ c„,^)^(F, = l) - < 5) 



where da 7T7T ' y (F n = 1) is the NLO cross section for e + e~ — > 7r + 7r~7 (initial 
state radiation only), inclusive in 8 n7T and 9 n under the assumption of pointlike 
pions, and corresponds to the quantity H of Eq. (2); A(s) is the fraction of 
7r + 7r _ 7i SR (7 FS R) events selected by the angular cuts 8 W7T < 15° or 9 n7T > 165°, 
50° < 6-n < 130° as a function of the invariant mass s of the virtual photon; 
and c S7vyS is a correction which must be applied due to the fact that, in the 
presence of simultaneous emission of initial- and final-state photons, s w is not 
identical to s, as it is in the case of ISR only. Both A(s) and c S7rjS have been 
obtained using the PHOKHARA Monte Carlo generator [13], which simulates 
the simultaneous emission of initial- and final-state photons. 

Note that <r(e + e~ — > 7r + 7r~) is obtained under the assumptions of (i) radiation 
emission from pointlike pions (the scalar QED model for FSR) and (ii) fac- 
torization, i.e., the absence of interference effects between the initial and final 
states [13]. We have used an alternative method which provides some test of 
the validity of the factorization ansatz and a valuable cross-check of the entire 
analysis. In this method, we correct the observed nrrj cross section for the 
relative amount of FSR expected from PHOKHARA, obtaining, in this way, 
a cross section that corresponds to ISR emission only. Next, we perform the 
event analysis, in which the acceptance correction and track-mass efficiency 
are taken from a Monte Carlo sample in which only ISR events are simulated. 
After dividing by the radiator function H, the full (i.e., real and virtual) FSR 
corrections to the cross section e + e" — > it + it~ are applied [26,27]. 
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The results for a(e + e~ — > 7r + 7r~) obtained with the two methods agree to 
within pa 0.2%. Taking into account the additional uncertainty arising from 
the assumption of radiation from pointlike pions, we assign an error of 0.3% 
due to the FSR corrections discussed in this section. 



3.2 Vacuum polarization corrections 



To obtain the pion form factor and the bare cross section, leptonic and hadronic 
vacuum polarization contributions in the photon propagator must be sub- 
tracted. This can be done by correcting the cross section for the running of a 
as follows: 

fa(0)\ 2 

Obare — CT dressed — — . (6) 



While the leptonic contribution Aai ep (s) can be analytically calculated, for 
the hadronic contribution, Aa had (s), we have used <Jh a d(s) values measured 
previously [28]. 

The pion form factor (F^s)! 2 obtained after additional subtraction of FSR 
is shown in Fig. 7. Note that in this case, since the FSR effects have been 
removed, s w = s. 

45 
40 
35 
30 
25 
20 
15 
10 
5 

0.4 0.5 0.6 0.7 0.8 0.9 

\ [GeV 2 ] 

Fig. 7. Pion form factor. 
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4 Results 



Our results are summarized in Table 1, which lists: 

• the differential cross section da(e + e~ — > 7i + 7i^'~f)/ds 7T as a function of the 
invariant mass of the di-pion system, s n , in the angular region 9 n7T < 15° or 
9 nn > 165°, 0° < 9 n < 180°; 

• the physical cross section a(e + e~ — > it + it~), which includes FSR and vac- 
uum polarization effects, as a function of the invariant mass of the virtual 
photon s; 

• the pion form factor with FSR and vacuum polarization effects removed, as 
a function of s (equal to s n in this case) . 

The errors given in Table 1 are statistical only, while the common systematic 
error is shown in Tables 2, 3, and 4. It should be noted that the statistical 
errors account only for the diagonal elements of the covariance matrix. The 
bin-by-bin errors are correlated as a result of the unfolding procedure; for error 
propagation, as for example in the calculation of a M (see below) , the covariance 
matrix must be used. 

The unfolding procedure is necessary in order to provide a table of data 
points at meaningful values of s n . However, the procedure itself introduces 
additional systematic uncertainties because of the numerical instability of the 
problem. For the comparison of our data with a specific theoretical prediction, 
we strongly recomend fitting our observed spectrum with a convolution of the 
theoretical curve and the detector response matrix, which is available upon 
request [29]. 

The <r(e + e~ — > 7r + 7r~) cross section, divided by the vacuum polarization, has 
been used to evaluate the contribution to d due to the tx + h" channel in the 
energy range 0.35 < s n < 0.95 GeV 2 . The resulting value (in 10~ 10 units) is 

0^(0.35, 0.95) = 388.7 ± 0.8 stat ± 3.5 syst ± 3.5 th . (7) 
The various contributions to the systematic error on are listed in Table 5. 



5 Conclusions 

We have measured the cross section for the process e + e~^ 7r + 7r~7 with the 
pion system emitted at small polar angles with respect to the electron or 
positron beam (9 nn < 15°, 9 n7T > 165°) in the energy region 0.35 < < 
0.95 GeV 2 . Using Eq. (2), we have derived the cross section for the process 
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e + e~^ 7r + 7r~, as listed in Table 1. These values, corrected for the vacuum 
polarization, can be used to derive part of the hadronic contribution to the 
muon anomalous magnetic moment with a negligible statistical error and with 
a systematic error of 0.9%(exp) © 0.9%(th). 



GeV 2 


nb/GeV 2 


7T + 7T 

nb 


\FJs)\ 2 


s n 
GeV 2 


7T + 7T 7 

nb/GeV 2 


7T + 7T 

nb 


\FJs)\ 2 





35 


13.40±0.24 


330±7 


7.68±0.16 





65 


59 


40±0 


28 


714±4 


24.69±0.15 





36 


14.59±0.24 


349±7 


8.26±0.16 





66 


56 


38±0 


24 


657±4 


23.05±0.14 





37 


15.78±0.24 


370±7 


8.92±0.16 





67 


53 


04±0 


23 


595±4 


21.18±0.13 





38 


17.04±0.24 


392±6 


9.60±0.16 





68 


49 


87±0 


26 


543±4 


19.57±0.13 





39 


18.63±0.23 


416±6 


10.35±0.15 





69 


46 


98±0 


22 


493.2±3.1 


18.02±0.11 





40 


20.34±0.27 


450±7 


11.40±0.17 





70 


44 


16±0 


21 


447.0±2.9 


16.54±0.11 





41 


22.64±0.24 


489±6 


12.59±0.16 





71 


41 


54±0 


19 


405.0±2.6 


15.17±0.10 





42 


24.56±0.27 


521±7 


13.63±0.18 





72 


39 


05±0 


21 


367.1±2.6 


13.92±0.10 





43 


27.07±0.28 


564±7 


15.01±0.18 





73 


36 


87±0 


17 


333.3±2.2 


12.78±0.08 





44 


29.99±0.27 


608±7 


16.43±0.18 





74 


35 


20±0 


18 


304.6±2.0 


11.81±0.08 





45 


32.65±0.28 


649±7 


17.82±0.19 





75 


33 


22±0 


16 


277.8±1.8 


10.89±0.07 





46 


36.24±0.27 


710±7 


19.79±0.18 





76 


31 


99±0 


16 


257.2±1.7 


10.19±0.07 





47 


40.10±0.29 


769±7 


21.78±0.20 





77 


30 


51±0 


17 


233.8±1.7 


9.37±0.07 





48 


44.34±0.31 


830±7 


23.86±0.20 





78 


29 


60±0 


16 


217.7±1.6 


8.82±0.06 





49 


48.94±0.28 


895±7 


26.11±0.20 





79 


28 


52±0 


13 


200.3±1.3 


8.20±0.05 





50 


54.1±0.4 


967±8 


28.60±0.23 





80 


27 


53±0 


14 


184.5±1.3 


7.63±0.05 





51 


59.77±0.32 


1041±7 


31.23±0.22 





81 


27 


00±0 


14 


172.1±1.2 


7.20±0.05 





52 


64.93±0.32 


1102±7 


33.50±0.22 





82 


26 


48±0 


13 


160.0±1.1 


6.76±0.05 





53 


70.24±0.35 


1171±8 


36.05±0.23 





83 


25 


84±0 


15 


148.4±1.1 


6.33±0.05 





54 


75.6±0.4 


1226±8 


38.20±0.26 





84 


25 


45±0 


13 


138.5±1.0 


5.97±0.04 





55 


80.2±0.4 


1279±8 


40.32±0.24 





85 


25 


16±0 


13 


129.2±0.9 


5.63±0.04 





56 


83.47±0.35 


1288±7 


41.07±0.24 





86 


24 


96±0 


12 


120.3±0.8 


5.29±0.04 





57 


86.06±0.34 


1302±7 


41.98±0.23 





87 


24 


81±0 


15 


111.8±0.9 


4.97±0.04 





58 


87.85±0.34 


1297±7 


42.36±0.23 





88 


25 


09±0 


14 


106.3±0.8 


4.774±0.035 





59 


89.5±0.4 


1282±7 


42.46±0.24 





89 


25 


17±0 


12 


99.5±0.7 


4.516±0.030 





60 


90.31±0.35 


1266±7 


42.58±0.23 





90 


25 


37±0 


13 


93.1±0.6 


4.269±0.030 





61 


74.20±0.35 


1006±6 


32.43±0.20 





91 


25 


86±0 


12 


87.6±0.6 


4.059±0.027 





62 


65.49±0.28 


857±5 


27.99±0.16 





92 


26 


87±0 


14 


83.0±0.6 


3.886±0.026 





63 


64.14±0.28 


817±5 


27.32±0.16 





93 


27 


94±0 


14 


79.1±0.5 


3.741±0.025 





64 


62.09±0.26 


772±4 


26.27±0.15 





94 


29 


49±0 


16 


75.3±0.5 


3.599±0.025 



Table 1 



Cross sections da(e + e~ — ► 7r + 7r _ 7)/(is 7r , a{e + e~ — > 7r + 7r~) and the pion form factor 
in 0.01 GeV 2 intervals where the value given indicates the lower bound. Note that 
while the 7r + 7r~7 cross section is given as a function of s^, the inr cross section and 
| 1 2 are given as functions of the invariant mass s of the intermediate photon 7*. 
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s (GeV 2 ) 





1 


2 


3 


4 


5 


6 


7 


8 


9 


0.3... 












0.8 


0.7 


0.6 


0.6 


0.5 


0.4... 


0.5 


0.4 


0.4 


0.4 


0.4 


0.4 


0.4 


0.4 


0.4 


0.4 


0.5... 


0.3 


0.3 


0.3 


0.3 


0.3 


0.3 


0.3 


0.3 


0.3 


0.3 


0.6... 


0.3 


0.3 


0.3 


0.3 


0.3 


0.3 


0.3 


0.3 


0.3 


0.3 


0.7... 


0.3 


0.2 


0.3 


0.3 


0.3 


0.3 


0.3 


0.2 


0.2 


0.2 


0.8... 


0.2 


0.2 


0.2 


0.2 


0.2 


0.2 


0.2 


0.2 


0.2 


0.2 


0.9... 


0.2 


0.2 


0.2 


0.2 


0.2 













Table 2 " " " 
Bin-by-bin correlated systematic error in % due to background subtraction in 0.01 
GeV 2 intervals. The indicated values for s represent the lower bin edge. 



s (GeV 2 ) 


0.58 


0.59 


0.6 


0.61 


0.62 


0.63 


0.64 


0.65 


&unf 


0.4 


0.9 


1.4 


3.6 


0.9 


0.8 


0.5 


0.4 



Table 3 ' 

Bin-by-bin correlated systematic error in % on da(e + e~ — > 7r + 7r~7)/(is 7r and 
cr(e + e~ — > 7r + 7r~) due to unfolding in 0.01 GeV 2 intervals. The indicated values 
for s represent the lower bin edge. 









IF I 2 


Acceptance 




0.3 % flat in s w 


Trigger 


exp(0.43 - 4.9s w [GeV 2 ]) % + 0.08 % 


Reconstruction Filter 




0.6 % flat in S7T 


Tracking 




0.3 % flat in SlT 


Vertex 




0.3 % flat in s n 


Particle ID 




0.1 % flat in S7T 


Trackmass 




0.2 % flat in s w 


Luminosity 
FSR resummation 
Radiation function (H(s n )) 
Vacuum Polarization 


0.6 % flat in s w 




0.3 % 




0.5 % 






0.2 % 



Table 4 

List of completely bin- by-bin correlated systematic effects. 
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Acceptance 


n q 07 
0.3 70 


Trigger 


0.3 70 


Reconstruction Filter 


0.O /o 


Irackmg 


0.3 /o 


Vertex 


0.3 70 


rarticle ID 


0.1 7o 


1 von /m occ 1 


U.Z /0 


Background subtraction 


0.3 % 


T Infoldine - 


0.2 % 


Total PYT) SVStPTTiatlCR 


0.9 % 


Luminosity 


n ft 0/ 
0.D 70 


Vacuum Polarization 


U.Z /o 


FSR resummation 


0.3 % 


Radiation function (H(s n )) 


0.5 % 


Total theory systematics 


0.9 % 



Table 5 

List of systematic errors on a M . 

Future improvements are expected using data taken in 2002, where more sta- 
ble background conditions and an improved trigger logic should allow for a 
considerable reduction of the systematic effects stemming from the offline re- 
construction filter and the trigger. A similar analysis, applied to events with 
9 W7r at larger angles, can probe the energy region down to threshold. Moreover, 
improved Monte Carlo generators both for the luminosity measurement and 
for the ISR process are expected to be available in the near future, which will 
help to reduce the theoretical contribution to the systematic error. 
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